Collective Quantum Dot Inversion and Amplification of Photon and Phonon Waves 
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The possibility of steady-state population inversion in a small sample of strongly driven two-level 
emitters like quantum dots (QDs) in micro-cavities, and its utilization towards amplification of light 
and acoustic waves is investigated theoretically. We find that inversion and absorption spectrum of 
photons, and phonons crucially depend on the interplay between the intrinsic vacuum and phonon 
environments. The absorption profiles of photons and phonons show marked novel features like 
gain instead of transparency and absorption reversed to gain, respectively. Furthermore, we report 
collectivity induced substantial enhancement of inversion and pronounced gain in the photon, and 
phonon absorption spectrum for a wavelength size QD ensemble. 

PACS numbers: 78.67.Hc, 43.35.Gk, 42.50.Ct, 42.50.Nn 



Even after fifty years of the advent of laser, achieving 
population inversion and lasing in novel systems remain 
a topic of continuing interest In particular, inversion 
in solid-state emitters like quantum dots (QDs) has at- 
tracted considerable interests recently @t3' due to the 
possibility of lasing in a driven two-level system con- 
trary to the common notion of multilevel manifolds 
Semiconductor QDs are novel nano-structures Q with 
unique engineerable features like large optical dipole mo- 
ments and emission wavelengths. They can also be in- 
tegrated efficient to micro-cavities and waveguides 0. 
Their similarities to two-level atoms have been demon- 
strated in phenomena like Autler-Townes doublet |9|, 
Mollow triplet in resonance fluorescence fil l [Til , spectral 
line narrowing [l2l - il4| and superradiance [15j]. Further- 
more, collectivity induced fluorescence and decoherence 
of a few nearly identical QDs have also been investigated 
16]. Additonally, due to potential application as solid 
state qubits numerous coherent optical studies of this 
system has been undertaken in the past decade [l7l - |2lT |. 
During last few years, sensitivity of QDs to the nature 
of their phonon environment has been exploited in nu- 
merous investigations. Phonon-assisted damping of Rabi 
oscillations [22h24|. e ffect of phonons on polarization- 
entangled photons [25] , phonon induced transitions of ex- 
citons to cavity photons [26] , heat pumping with optically 
driven excitons in phonon environment [271 ] and effect of 
electr on-p honon coupling on resonance fluorescence spec- 
trum [28-30] have been reported. Furthermore, phonon 
mediated population inversion in a single QD-cavity sys- 
tem driven by a continuous wave [3l| and a pulse laser 
excitation [321 ] was proposed recently. 

In light of these works, here, we propose a novel scheme to 
achieve and enhance population inversion in an ensemble 
of strongly driven two-level solid-state emitters like QDs 
via coupling to a phonon-photon reservoir. We find im- 



portant interplay of the phonon-photon reservoirs on the 
QD-dynamics and in creation of steady-state inversion in 
both bare- and dressed-states. As a key finding, we re- 
port phonon induced unique features in the QD photon 
absorption spectrum like, gain instead of transparency 
at resonance and absorption switched to gain for bare- 
state inversion. We also find gain in the phonon absorp- 
tion spectrum under condition of dressed-state inversion 
which can be harnessed for amplification of phonon waves 
at THz frequencies. Additionally, for an wavelength size 
ensemble of QDs, we report collectivity assisted enhance- 
ment of inversion and improvement of gain in the photon 
or phonon absorption spectra. 

Model: We consider an ensemble of QDs embedded in a 
substrate (such as InGaAs encased in GaAs substrate) or 
in a micro-cavity (like that formed by InAs/GaAs Bragg 
reflectors) illuminated by a laser of frequency w^. As 
shown schematically in Fig. ([I} the ground state |1) and 
the state containing a trapped electron-hole pair (exci- 
ton) 1 2) in each QD form a two- level system with transi- 
tion frequency uj x . The linear dimension of the ensemble 
is smaller or of the order of the relevant emission wave- 
length. The laser is detuned from the exciton transition 
by A = uj x — ujl and drives it with a Rabi frequency of 
il. Furthermore, the QDs interact with the solid state 
environment (phonon reservoir) and with field modes of 
the electromagnetic vacuum (photon reservoir) leading 
to different incoherent decay processes. 
The Hamiltonian describing the interaction of N QDs 
with the laser and the vacuum and phonon reservoirs 
in a frame rotating with the laser frequency wj, can be 
represented as: 

N 
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FIG. 1: (color online) Schematics of interaction of a collection 
of QDs with a coherent field and phonon reservoir, (a) QDs 
as two-level systems embedded in a substrate or in a micro- 
cavity of emission wavelength size, driven by a laser of Rabi 
frequency f2. (b) In strong field, the QDs are dressed by the 
laser giving rise to dressed-states |±) and characteristic Mol- 
low triplet in emitted photon radiation. However in presence 
of phonons, additional transitions at the Rabi frequency 2f2 
among the dressed-states are induced. 



In Eq. ([T]) the first three terms represent the free Hamil- 
tonians of the electromagnetic vacuum, the phonon reser- 
voir and the two-level QDs. The fourth term corresponds 
to the laser-QD interaction, while the two terms inside 
the curly brackets describe, respectively, the QDs inter- 
action with the surrounding phonon and vacuum reser- 
voirs. Both the vacuum and phonon reservoirs are rep- 
resented in the form of harmonic oscillators with a k (a,k) 
and b p (b p ) as the creation (annihilation) operator of the 

vacuum held in mode k and phonons in mode p hav- 
ing frequencies and lu p respectively. They satisfy the 
standard boson commutation relations: [Q u ,Q u i] = fiuu' 
and [Q U ,Q U .] = [QtQl] = 0, (Q = {a,b},u= {k,p}). 
The laser is detuned to the vacuum mode k by Sk = 
uik — ojl. The terms A p , cfk and pj are respectively the 
QD- phonon coupling constant, the mode function of the 
three-dimensional multi-mode field (33j | and the dipole 
moment of the j-th QD. The raising and lowering op- 
erators for each QD are denoted by Sj = \2)jj(l\ and 
Sj = (S+)t, while S| = (\2) n (2\ - |l)„(l|)/2, obey- 
ing the standard su(2) angular momentum commutation 
relations. 

In the following we consider the QDs system driven by 
an intense laser. This leads to the dressing of both 
the bare-states in each QD creating the dressed-states, 
= sin0|l)j + cos0|2)j, = cos6\l)j - sin9\2)j 
where tan 29 = (2J7/A). Hence, it is imperative to con- 
tinue the further analysis in the dressed-state basis. To 
study the collective quantum dynamics of the strongly 
driven QD ensemble interacting with vacuum and phonon 
reservoirs, we next derive a dressed basis master equation 
(DME) in the interaction picture. 

Collective Quantum Dynamics of QDs: The general 



form of the master equation in the interaction pic- 
ture is given by [Hj]: dp sb /dt = l/ih[Hi(t), p sb (0)] - 

l/h 2 Sodt'[Hi(t),[Hi(f),Psb{t')]]. Here, p sb is the den- 
sity operator of the combined QDs + reservoirs (both 
vacuum and phonons) in the interaction picture. The 
dressed-state interaction Hamiltonian Hi can be evalu- 
ated by using the dressed-state transformations in Eq ([1]) 
and an unitary transformation by the operator e ( l /?>-) H ot j 

where H = Y,k h5kd\a k + Y, P ^P b l b P + ^X)jli R j- 
On substituting the interaction Hamiltonian and taking 
trace over the vacuum field and phonon modes, after 
some tedious algebra, we get the DME for the reduced 
density operator p of the QDs in the Born-Markov ap- 
proximation, and weak coupling regime as: 

^ + in[R*, P } = -xoJ2[ R i> R lp\ 

1,3 

-x+ £ [Rt,R]p] -x-J2 [ R T' R tp\ + H - c -P) 

1,3 1,3 

where VL = fi - A p with Q = ^(A/2) 2 + fl 2 be- 
ing the generalized Rabi frequency, while A p is an ad- 
ditional shift due to QD-phonon coupling. The well 
known Lamb shift due to vacuum fields has been ab- 
sorbed in the exciton frequency. The dressed basis op- 
erators are defined as: i?| = — |— | and 
Rj = \+)jj{—\,Rj — (R^f V obeying the commutator 
relation [B^ p> R%,] = S jt {6 ^R% - Sp, a R%}. 
In Eq. ((2|) xo,± are the renormalized decay rates 
given by: X o = (T sin 2 20 + T d cos 2 26) /4, X + = 
[(n+ 1)T P + r d ]sin 2 26>/4 + r+ cos 4 9 and 
X- = (nr p + r d ) sin 2 29/4 + r_sin 4 6». Here, 

r = (g/h) 2 K/[ K 2 + s 2 c ], r± - (g/h) 2 K /{K 2 + (5 cT 2n)% 

and T p = r p (2ft) = tt J2 p {K/ h f %p - 2n )- Thc 
term Q = g ■ p is the vacuum-QD coupling constant 
while sin20_ = fi/Q, cos 2 9 = (fi + A/2)/2f\ and 
sin 2 6» = (fi - A/2)/2fi. Further, n = n{2Q) = 
[exp (2hVL/k B T) - I]- 1 is the mean phonon number at 
frequency 2S7 and temperature T, with ks being the 
Bolzmann constant. The simplified expressions for Fo ; ± 
is derived under the assumption that the QDs interact 
with a broad-band cavity mode of frequency lo c and 
decay rate k. The driving laser is detuned to the cavity 
frequency by S c = ui c — uj^. Note that the correlation 
functions for the cavity vacuum field and phonon 
operators are taken in the standard form [33]. Thc 
decay terms included in the renormalized decay rates 
can be physically explained as follows: Tq and T± are 
the radiative decay rates of the QDs dressed transitions 
involving the upper and lower doublets |±) — > |±) and 
|±) — > |=f) (|2) — > |1)), respectively. T p is the phonon 
induced dissipation rate of the QDs at the generalized 
Rabi frequency 2ti between the dressed-state doublets 
of |2) as well as |1) (see Fig. lb) and Td is the pure 
dephasing rate of QDs due to phenomenon like elastic 
phonon scattering [25| . Note that in deriving Eq. @ 
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in the standard Lindbald form, we have performed the 
secular approximation, i.e. we neglected the terms 
rotating with the frequency 2ft and higher meaning that 
2ft ^> AT(1 + n). Here T _1 typically corresponds to the 
smallest time scale of the system. The corrections to 
the results obtained in the secular approximation are of 
the order of [iVT(l + n)/ft] 2 and can be neglected in the 
intense field limit. 

Fluorescence and absorption spectrum: To understand 
the underlying physics of the QD dynamics, we next 
investigate the steady-state fluorescence and absorption 
spectrum of a single QD. The total steady-state photon 
fluorescence spectrum defined in the standard form (33[, 
under secular approximation and in the dressed-state ba- 
sis is evaluated as: St{uj) = S c (u>) + S(oj), where S c (ui) 
and S(u) corresponds to the coherent and incoherent 
parts and are given by: 



S c (u) = -{R Z ) 2 SS sin 2 2^(w- WL ), 



S(u) = sin 4 
+ sin 2 29 

+ cos 4 e 



(3) 



ML 



X++X-J V 2 ML + {u-u L + 2ft) 2 



ML 

Fmc 



(X+ + X-) 2 T 2 MC + (u: - u L y 

X- \ r A/i?, 



X++X-J T 2 MR + 2ft) 2 



(4) 



The corresponding dressed-state absorption spectrum 
[34| in the secular approximation, for photons is : 



A(v) = i(R z ) ss { sin 



} Tml + i(y_ - u L + 2ft) 
T 2 ML + {v-u L + 2h) 2 

4 Tmr + ~ u L - 2ft) 1 

cos 9 = — >, (5) 

T 2 MR + {v-u JL -2tl) 2 y 



while for phonons is given by: 



A(y ) = -L(IV)„ sin 2 2g( ^ + i[Vp r 

P 4 lr 2 /i + K~2ft) 2 l 



2ft) 1 



(6) 



Here, (R z ) ss is the single QD dressed basis steady-state 
population inversion given by: — x+)/(X- + X+)- 
Ta/l = T RL = 4xo + X+ + X- and ^mc = 2(x+ + 
X—) with is, Up being the weak photon and phonon probe 
frequencies, respectively. In Fig. [3Ja) & [^b) we plot 
Eq. (|4]) and imaginary part of Eq. ([5]). The solid and 
dashed lines in the plot corresponds to resonant and non- 
resonant coherent excitation of the QD. We have assumed 
bad cavity limit i.e. k>5,A p and approximated ro,± = 
7, with 7 being the QD radiative decay rate in the cavity. 
The photon fluorescence spectrum shown in Fig. HJa) is 
similar to the Mollow spectrum well known for atoms, 
however, with certain key modifications due to the 
phonon bath. The position of the Mollow sidebands are 
shifted to ±2ft with phonon induced broadening of their 
widths Tml and that of the central line Tmc as can be 



easily seen from their respective analytical forms. Addi- 
tionally, in presence of phonons the sidebands are asym- 
metric with the ratio of the right to left peak heights as: 
[(2ft + A)/(2ft - A)] 2 x-/x+- This asymmetry is pro- 
nounced for non-resonant driving of the QD as shown in 
Fig. HJa) by the dashed curve. For resonant drive how- 
ever, the asymmetry can be prominent in the conditions: 
n <C 1 (low temperature reservoir ) and T p ;§> (1^,7). 
The coherent part of the fluorescence spectrum S c (uj), 
that signifies the elastic scattering process in presence of 
phonons, does not vanish even at resonance and has an 
amplitude oc r 2 /[(2n + l)r p + 2(r d + 7)] 2 . Note that, 
such phonon induced modifications of the Mollow triplet 
has been observed in a recent experiment [30| . 
The imaginary part of the photon absorption spectrum 
shown in Fig. [U[b) is strikingly different from a typical 
Mollow absorption spectrum of a two-level atom. As a 
key finding, for resonant excitation of a QD we report 
phonon induced gain and absorption in the spectrum at 
v — lol = =F2ft respectively with peak heights oc {R z ) ss /A. 
Note that, these characteristics are absent in atomic sys- 
tem, where there is no phonon coupling and the probe 
field is transparent at resonance as (R z ) ss = 0. Further- 
more, for both resonant and non-resonant excitation we 
find from Eq. (J5j> broadening of the gain and absorption 
lines due to phonons. The dispersion characteristics of 
the absorption spectrum given by the real part of Eq. ([5]) 
behave in accordance with the absorption characteristics 
and have no additional features. 

The imaginary part of Eq. © exhibits another remark- 
able feature, gain in the phonon absorption spectrum for 
population inversion in the dressed-state. This gain can, 
in principle, lead to amplification of a weak phonon wave 
in presence of an additional phonon cavity [3ol I3H ] tuned 
to the dressed transition at the shifted generalized Rabi 
frequency of 2ft (see Fig. lb). Note that the phonon 
gain profile is enhanced only in presence of the vacuum 
reservoir as can be seen from the expressions for xo,±- 
For a collection of N QDs in an ensemble of the size of 
emission wavelength, the emission-absorption spectrum 
is modified due to collective effects. The spectral lines 
are broadened and enhanced, with the enhancement pro- 
portional to A^ 2 and N, respectively, for the Mollow peaks 
in photon fluorescence and gain profile in the photon or 
phonon absorption [34]. In the following we discuss an- 
other key finding of the paper: collectivity assisted en- 
hancement of steady-state population inversion. 
Population inversion in QDs: In order to study the col- 
lective steady-state population inversion of the QDs in 
bare and dressed basis we need to obtain the steady- 
state solution of Eq. ([2]). For this purpose we first de- 
fine the dressed-state collective operators as: IZ^ = 
Y^j=i^f^ z = Ylj=iRj and then look for a solution 
of the form tol 



z- 1 exp[-nK z ], 



(7) 



where the normalization Z is determined by the require- 
ment Tr{/9 S } = 1. The unknown variable r\ can be eval- 
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FIG. 2: (color online) The fluorescence and absorption spec- 
trum of a QD at resonance (solid curve) and for negative 
detuning (dashed curve) of A = —40 /ieV. The y-axis in both 
the plots are in dimensionless units and in (b) is multiplied 
by a factor of 10. We have considered parameters similar to 
that achievable in experiments: £7 = 80 ueV, A p = 20 peV, 
r = r± = 7 = 0.84 fieV, T p = 0.34 fieV, F d = 3.84 ueV and 
T = 6K. 

uated by inserting Eq. ([7]) in the DME and taking the 
steady-state dp/dt = condition. On doing so we ob- 
tain, 

V = ln( X+ /x-)/2. (8) 

We next look for the inversion in the bare-states. In the 
secular approximation the collective bare-state inversion 
is given by, 

(S z ) ss = (A/2fi) (R z ) ss /2, (9) 

where (1Z Z ) SS is the collective steady-state inversion in 
the dressed-state found from Eqs Q and (jSJ as @: 

/ l + eW \ 2 (eW 1 )-e 2 ") 
K >ss U - eWN+V J (1 - e 2 ")(l - e^^+D) ' 

To gain physical insight into the behavior of the bare- 
state population we need to understand the inversion of 
dressed-state populations. For this purpose we consider 
two limits of Eq. ©: (i) e 2 '' = f = X+/X- < 1 and ( n ) 
£ = x+/X- ^ 1- From Eq. (fT0|) . we get in the limit (i): 
(TZ Z ) SS /N ~ 1 - [2£/iV(l - £)], which clearly approaches 
maximum inversion (H z ) ss /N = 1, for JV 3> 1. As iV 
represent the number of QDs in the ensemble, we can 
conclude that collectivity among the QDs enhances pop- 
ulation inversion in the dressed-state. In the limit (ii) on 
the other hand we get: (TZ Z ) SS /N ~ [2/iV(£-l)]-l which 
for JV > 1 becomes, (1Z Z ) SS /N = —1, suggesting that 
all the population tends to reside in the lower dressed- 
state |— ) and, thus, no inversion. Note that, in the limit 
X+/x- 1? (7l z ) ss /N — > there is also no inversion 
as the dressed-states are equally populated. Addition- 
ally, we would like to emphasize that inversion is possi- 
ble only in presence of the radiative decay T± among the 
dressed-states with T± > T p . This can be proved by an- 
alyzing the expressions for \±- On substituting T± = 
wehavex+/x- = [{n+l)T p +T d }/[nT p + r d } > 1, for any 
bath temperature and QD parameters thereby implying 
no possible dressed-state inversion. 

It is clear from the above discussion and Eqs © and 
^ , that for conditions of inversion in the dressed-state 
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FIG. 3: (color online) (a) Steady-state population inversion 
in the bare-state for N — 10 (dotted curve), N = 5 (dashed 
line) and N — 1 (solid curve) QDs. (b) Absorption spec- 
trum of a single QD for A/2Q, — —1 (maximum inversion 
in Fig. 3a) in presence (solid curve) and absence (dot-dashed 
line) of phonons. In (b) the y-axis is in dimensionless units 
and multiplied by a factor of 100 while the x-axis by 0.01. 
Parameters for the plots are: Q — 2 meV, A p = 20 /xeV, 
r = r± = 7 = 0.84 ueV, F p = 207, r d = 407 and T = 6K. 

there will be no population inversion among the bare- 
states |1) and |2). However, at low temperature (small 
rl), if T p ^> r± we can satisfy the condition x+/X- ^ 1 
and (1Z Z ) /N < thereby leading to inversion among the 
bare-states for negative laser detuning. 
In Fig. IS^a), we show the steady-state population inver- 
sion among the bare-states for different number of QDs 
and typical experimental parameters [13] ■ Inversion is 
seen to be achievable for certain range of negative laser 
detuning as expected from the above discussion. Ad- 
ditionally, we find collectivity induced enhancement of 
population inversion in the bare-state. This follows di- 
rectly from the behavior of {1Z Z ) SS /N for x+/X- ^ 1- 
We report almost 40% inversion for reasonable (JV = 10) 
number of QDs (complete inversion in bare-state corre- 
sponds to a (S z )/N = 0.5). Note that the mean-phonon 
number depends on the generalized Rabi frequency and, 
thus, varies as detuning varies. This has been explicitly 
considered for the plots and we have checked that the 
dependence does not affect our findings and conclusions 
for the domain of parameters used. 

In Fig. EJb) we show the behavior of the imaginary part 
of photon absorption spectrum of Eq. (J5j) at maximum 
bare-state inversion for a single QD. The spectrum show 
completely new feature in presence of phonons - gain at 
bare-state inversion. However, for r± S> T p or in absence 
of phonons, we get back the typical absorption peak in- 
stead of gain (dot-dashed line in Fig. [2b). This implies 
that we can get amplification of a probe laser from a QD 
ensemble. 

Applications: For a small ensemble of QDs in a micro- 
cavity, we find that collectivity enhances the inversion 
substantially both in bare and dressed-state. Hence, the 
phonon induced gain characteristics of photon-absorption 
spectrum that we report here for bare-state inversion, can 
be utilized towards amplification of a weak light field at 
optical frequencies in such solid-state system. The pho- 
ton induced gain characteristic of the phonon-absorption 
spectrum for dressed-state inversion, on the other hand, 
can amplify a phonon wave at THz frequencies tuned to 
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the 2f2 transition (as ~ ps is achievable in QDs). 
In conclusion, we have derived a collective dressed-state 
master equation for an ensemble of strongly driven QDs 
coupled to phonon and vacuum reservoirs. As a key find- 
ing, we reported gain characteristic in the photon and 
phonon absorption spectrum instead of transparency and 
absorption respectively due to crucial interplay of the 
vacuum and phonon reservoirs. Furthermore, we have 
analyzed the conditions for achieving population inver- 
sion in the dressed and bare-states. We have shown that 
collectivity among the QDs can enhance the obtained 



population inversion. Finally, as applications, we have 
proposed amplification of phonon waves at THz frequen- 
cies and light fields at optical frequencies. Our find- 
ings show the possibility for amplification of photon and 
phonon waves in QDs and can open new directions to- 
wards creation of nanoscale optical and acoustic coherent 
sources. 
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